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Abstract 

In wireless relay networks, noise at the relays can be correlated possibly due to common interference or noise 
propagation from preceding hops. In this work we consider a parallel relay network with noise correlation. For the 
relay strategy of amplify-and-forward (AF), we determine the optimal rate maximizing relay gains when correlation 
knowledge is available at the relays. The effect of correlation on the performance of the relay networks is analyzed 
for the cases where full knowledge of correlation is available at the relays and when there is no knowledge about the 
correlation structure. Interestingly we find that, on the average, noise correlation is beneficial regardless of whether 
the relays know the noise covariance matrix or not. However, the knowledge of correlation can greatly improve the 
performance. Typically, the performance improvement from correlation knowledge increases with the relay power 
and the number of relays. With perfect correlation knowledge the system is capable of canceling interference if the 
number of interferers is less than the number of relays. 

For a dual-hop multiple access parallel network, we obtain closed form expressions for the maximum sum-rate 
and the optimal relay strategy. The relay optimization for networks with three hops is also considered. For any relay 
gains for the first stage relays, this represents a parallel relay network with correlated noise. Based on the result of 
two hop networks with noise correlation, we propose an algorithm for solving the relay optimization problem for 
three-hop networks. 
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I. Introduction 



Wireless mesh networks where information is transferred via multiple hops and routes provide significant 
throughput enhancement and have been the focus of much recent work [l]-[4]. Various relay strategies have 
been studied in literature [5]-[7]. These strategies include amplify-and-forward [8], [9], where the relay sends a 
scaled version of its received signal to the destination, demodulate-and-forward [9] in which the relay demodulates 
individual symbols and retransmits, decode-and-forward [5], [7] in which the relay decodes the entire message, 
re-encodes it and re-transmits it to the destination, and compress-and-forward [5], [6] where the relay sends a 
quantized version of its received signal. Among these strategies, amplify and forward has been found to be highly 
suitable for parallel relay networks for its ability to pass on soft information [10]-[12]. In one of the first works [13] 
on the parallel relay channel, it is shown that AF achieves the capacity cut-set bound at high relay power. In [4] it is 
shown that full degrees of freedom can be achieved with AF in a multi-hop parallel relay network where the hops 
are orthogonalized by design. There are numerous other examples in literature that corroborate the effectiveness of 
AF [12], [14]-[19]. For an AF relay network, the relay design involves optimizing the relay amplification factors 
to maximize performance. Previous work pertaining to relay optimization for AF relay networks includes both 
single-user relaying [20]-[25] and multi-user relaying [2], [26]. 

Most work in the literature [20]-[27] assumes independent noise at the relay terminals. However, noise correlation 
between nodes occurs in wireless relay networks due to several reasons. In this paper we explore the effect of noise 
correlation for the following two models. 




Effective noise at the relays 



vi = aiii + 0212 + ni 



V2 = fciii + 62*2 + "-2 



Fig. 1. Noise correlation model 1 : Parallel relay network with common interference 



1) Common Interference Model 

Due to the broadcast nature of wireless networks, the relays are exposed to a set of common interferers 
resulting in correlated noise at the nodes. As shown in Fig. [H each relay in addition to its local noise observes 
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common interference. 
2) Noise propagation with multi-hop AF relaying (Three-hop model) 

Consider a three-hop network as shown in Fig. [2l If the first stage relays amplify their received signal, each 
relay in the second stage observes a linear combination of noise from the preceding stage relays along with 
its own local noise. This clearly results in correlated noise at the second-stage relays. 




First stage relay gains: di, d2 

Effective noise at the second stage relays 

vi = dihiini + d2h2in2 + gi 
V2 = dihi2ni + a!2'i22"2 + <?2 



n2 92 

Fig. 2. Noise correlation model 2: Three-hop parallel relay network 

Both the above models are of considerable practical importance. The natural question is whether the relays can 
exploit the correlation structure to improve performance. In practice, learning correlation may result in network 
overheads. Whether such overheads are justified depends on the potential advantage of learning correlation. Thus, 
our goal in this work is to estimate the improvement in performance when perfect correlation knowledge is available 
at the relays. The optimal relay design in this case will have the following two objectives: 

1) Increase signal power: The relay gains can be designed such that the copies of desired signal adds up in phase. 

2) Reduce interference power: The relay gains may be chosen such that the common interference terms add out 
of phase. 

In general, these objectives cannot be achieved simultaneously and the relay design can be expected to be a trade-off 
between the two objectives. 

This paper also addresses the relay optimization problem for networks with three hops. Such a network is 
interesting as it combines multiple hops with multiple paths. The problem of routing in the traditional link based 
decode and forward involves selecting a relay at each hop that maximizes the end to end throughput between the 
original source and the final destination. If there are N relays in each of the M hops, the problem is to find the 
best route among N^^'^~^'> routes. With AF, all the routes are utilized with priority for each route determined by the 
power allocation and phase equalization at each of the relays across multiple hops. As we can control the priority 
for the routes, optimizing relay gains may be viewed as soft-routing. 
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A. Results 

The main results of this paper are summarized below. 

1) We obtain closed form solutions for the optimal relay amplification vector and the maximum sum rate for 
the network with correlated noise at the relays. This result generalizes the single-user relay optimization in 
[20], [21] and the multi-source relay optimization in [2] both of which assume independent and identically 
distributed (i.i.d.) noise at the relays. 

2) We find that correlation, on the average, is always beneficial regardless of the presence or absence of correlation 
knowledge at the relays. This is true irrespective of the channel state information (CSI) at the relays. 

3) We study the benefits of exploiting noise correlation at the relays. We compare the maximum sum rate without 
correlation knowledge (relays use optimal amplification factors based on uncorrelated noise assumption even 
if noise is correlated) versus capacity with correlation knowledge. The following key questions are answered: 
Does correlation help? Is the correlation knowledge more (or less) helpful as the number of relays increases? 
What is the effect of correlation as the power at the relays increases, and when the first hop/second hop 
becomes stronger? We also provide asymptotic results to characterize the impact of relay noise correlation. 

4) We apply the results of two hop relay networks with correlated noise to solve the relay optimization problem 
for three-hop AF relay networks with independent noise processes at the relays. We also characterize the 
behavior of the three-hop network as the transmit power at any of the stages tends to oo. When the first or 
the last stage power is very high, the network can be reduced to a two-hop channel. When both the first and 
last stage powers are high, the network can be reduced to a point to point MIMO channel with a single data 
stream. 

B. Notation 

We use bold upper letters to denote matrices and bold lower letters to denote vectors. Further (.)*, (.)-^, {.)^ stand 
for conjugation, transposition and Hermitian transposition, respectively. a0b represents element wise multiplication 
of two vectors. A = diag(a) denotes a diagonal matrix A that contains on its diagonals the elements of a. det(A), 
Tr(A), Aniax(A), and Vmax(A) denote the determinant, trace, principal eigen-value, and principal eigen- vector of 
A respectively. A^^" denotes the Hermitian transposition of A^-*^. v(i) denotes the i^^ element of the vector v 
while Ajj denotes the element in the i*^ row and j^^ column of the matrix A. E stands for the expectation operator. 
V ~ CAA(w, A) indicates v is a complex Gaussian column vector with E[v] = w and E[vv1^] = A. 

This paper is organized as follows. In the next section we consider a two hop relay network with noise correlation 
and address the relay optimization problems for both single and multi-user scenarios. We analyze the impact of 
correlation in Section HHI In Section |IV] we address the relay optimization problem for three-hop AF networks. We 



Fig. 3. Two-hop single user parallel relay network 



conclude with Section |Vl 

II. Two-hop parallel relay network with noise correlation 

A. System Model 

We consider a two-hop parallel relay network as shown in Fig. [3] In this model, the source S communicates to the 
destination D through a set of N half-duplex relays, Ri, i = I...N. The data transmission takes place in two time 
slots. In the first slot, the source transmits to the relays and in the second slot, all the relays simultaneously forward 
their received signal to the destination. Note that there is no direct link between the source and the destination. All 
the nodes are equipped with a single antenna. The x 1 channel between the source and the relays is f while 
g is the X 1 channel between the relays and the destination. The entries of f and g are independent and i.i.d. 
complex Gaussian random variables with zero mean and unit variance. The source has power P and the relays 
have a total power of Pr. The relay received symbols during the first time slot are 

r = fx + nR (1) 

where x is the source transmitted signal with power <f = P, and hr ~ CM{0, K) is AWGN with the 
covariance matrix K, given by 

K = E[nRnJ^]. 

In the second slot, the relay Ri}^i scales its received signal with a complex scaling factor di and transmits to the 
destination. We collect the relay gain factors in an A^ x 1 vector d. Then the received signal at the destination 
can be expressed as 

y = d^Gix + d'^GnR + no (2) 

where G = diag(g) and no CM{0, 1) is additive white Gaussian noise (AWGN). Further, we assume that there 
is no correlation between destination noise n£, and the relay noise vector riR as these noise processes occur at 



two different time slots. We denote this relay network with the shorthand notation 9^2(^', f , d, g) where the 
subscript 2 indicates that the network consists of two hops. The sum power constraint of the relays can be expressed 
as 

dt [ffftp + K) Oil d = PR (3) 



where represents the Hadamard product or the element-wise multiplication. Note that there is a strict equality 
for the power constraint. This is because higher total relay transmit power is equivalent to lower noise power 
at the destination. Finally a word about the sum-power constraint. With individual power constraints, not all the 
relays may operate at full power [28]. Thus the total power expenditure at the relays may change depending on 
the channel conditions. With sum-power constraint relays use all the power. This is helpful in interference limited 
systems where the interference from the relays to other wireless links in the network needs to be estimated. And 
most importantly, the sum power constraint allows tractable analysis and useful insights can be obtained as a result 
of it. 

B. Relay Optimization 

From the destination received signal in ([2]), the signal-to-noise-ratio (SNR) at the destination is given by 

and the corresponding transmission ratq^ is i? = log(l + SNR). We seek to maximize the transmission rate with 
respect to the amplification vector d subject to the sum power constraint at the relays. The following theorem 
provides the optimal relay design and maximum achievable rate with amplify and forward relays in a two hop 
parallel relay network with noise correlation. 

Theorem 1: The maximum achievable rate in the two-hop parallel AF relay network 9^2(^'5 f , d, Pr, g) where 
the relay noise covariance matrix is K, is given by 

i?° = log(l + SNR°) (5) 

SNR° = PPH(f 0g)^A-i(f0g) (6) 

and the optimal relay amplification vector d° is 

d° = K(A-i(f0g))* (7) 
'Gaussian inputs at the source are optimal for AF relays. 
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where 



A = K ggTPR + ff^P I + K I 



(8) 



and 



(9) 



(f g)tA-i [(fftP + K) I] A-t(f g) ■ 



The proof is presented in the Appendix. 

Remark: In the absence of correlation, A is diagonal and the phase of relay Ri is such that it cancels the phase 
of the forward channel g(i) and the backward channel f(z). However, with noise correlation, A is not a diagonal 
matrix. Therefore the optimal relay gain may not co-phase the input signal as part of noise can be canceled with 
relay phase and gain adjustments. 

C. Properties of the optimal relay network 

In this section, we consider the hmiting properties of the optimal relaying scheme. The proofs are available in 
the Appendix. 

I) Relay power Pr oo: When the relay power is high, the network is equivalent to a single-input-multiple- 
output (SIMO) system modeled by y = fx + n where K is the noise co variance matrix at the multiple antenna 
receiver and P is the transmit power of the single antenna transmitter. Its capacity is given by 



Note that the system is independent of g. Here d° = K(G^TK^^f)* achieves (fTOb . When K is not invertible 
(or equivalently det(K) = 0), the capacity is infinite. This means that the effective noise at the destination can 
be completely eliminated regardless of the source transmission rate. For any relay gain vector, d, the network 
is equivalent to a SIMO system with the receive combining vector given by (d g). 

2) Source power P ^ oo: For the optimal d°, the network is equivalent to a MISO system represented by 
y = gx + n where n is unit variance AWGN. Its capacity is given by 



which is independent of the relay noise covariance matrix K. d(i) = Kg*(i)/f(i) achieves ([TTI ) in the 
relay network. For a general d, the relay network represents a multiple-input-single-output (MISO) system 



C = log(l + f^K"^fP). 



(10) 



C7 = log(l+gtgP^,) 



(11) 



y = gx + n, with the beamforming vector 



(dOf) 



and transmit power equal to Pr. 



V(d0f)t(d0f) 
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fL 



Fig. 4. Two-hop multi-source parallel relay network 

3) Singular K or det(K) = 0: Singula]] K in any of the noise correlation models of common interference and 
three-hop network may be unrealistic as the relays have independent local noise components. However the 
case of singular K provides useful insights when the power of the interference is very high for the common 
interference model, and when the total power of the first stage relays is very high in the three-hop model. 
Notice that the rank of K indicates the number of independent noise sources at the relays. Thus for singular 
K at most (N — rank(K)) noise terms can be canceled with appropriate relay gains. Due to destination noise, 
the capacity of the relay network is finite even when K is singular. 
The relay optimization result of Theorem 1 can be extended to several multi-user scenarios. As an example we 
consider the multiple access relay channel in the following. 

D. Extension: Multiple access parallel relay network with noise correlation 

Consider a multi-source parallel relay network as shown in Fig. ID In this model, L source nodes wish to 
communicate to a common destination with the assistance of N AF relays. In the first slot all the sources transmit. 
The relay received signals are given by 

L 

r = ^fiXi + nR (12) 

i=l 

where Xi is the signal transmitted by the i source whose power is Ellxil"^] = Pi. As in the single user case, 
the relays scale their received signal and transmit to the destination in the second slot. The received signal at the 
destination can be expressed as 

L 

y = d^Gf.x, + d'^GnR + no (13) 

i=l 

The following theorem provides the maximum achievable sum rate and the optimal relay design. 

Theorem 2: The maximum achievable sum rate in a two-hop multi-source parallel AF relay network with noise 
^Theorem 1 holds for any noise covariance matrix and we do not assume any structure for the covariance matrix. 
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correlation is given by 



i?° = log(l + SNR°) 



(14) 



SNR° = PnXr^^ [A-iR] 



(15) 



and the optimal relay amplification vector is 




Kax [A-^B])* 



(16) 



where 



L 



A = K ggtp^ + ^f^ftp^ 1 + K0 I 



(17) 



k=l 



L 



B = 5^(fk0g)(fk0g)^Pfe 



(18) 



fc=i 



and K ensures compliance with the relay power constraint. 

Refer to the Appendix for the proof. Since our principal goal is to investigate the effect of correlation in a single 
user parallel relay network, the rest of the paper will focus only on the single-user case. 



In this section, we address the following important questions: 

1) Does correlation help when the relays know the co variance matrix? 

2) Does correlation hurt when the relays are unaware of it? 

To answer the first question we compare the following two scenarios. 

• The noise is correlated and relays are aware of the correlation. (Scheme- 11) 

• The noise is uncorrected. (Scheme-00) 

To answer the second question we compare the following two scenarios. 

• The noise is correlated and relays are unaware of the correlation. (Scheme- 10) 

• The noise is uncorrelated. (Scheme-00) 

If Scheme- 1 1 outperforms Scheme-00 then we can say that correlation helps if the relays are aware of it. Similarly, 
if Scheme-00 outperforms Scheme- 10 then we can say that correlation hurts when the relays are unaware of it. 



III. Does Correlation Help? 



A. Benchmark Schemes 

1) Relays with uncorrelated noise: Scheme-00: This setup has been commonly studied in literature [2], [20], 
[21], [24]. The system model consists of a two hop parallel relay network with independent noise at the relays. 
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The noise covariance matrix, which is diagonal, is given by K I. For a sum power constraint at the relays, the 
optimum relay amplification vector is found in [20], [21]. The result can also be obtained from Theorem 1 by 
replacing K with K I. We thus obtain 



doo = n 



(GKGtPfl + ff tp + K) ll (f^0g^f 



(19) 



The SNR achieved at the destination with the optimal relay gain doo is given by 



SNRoo = PPii(f g)^ [(GKGtp^ + fftp + K) ll (f g). 



(20) 



A special case of Scheme-00 is where the relay noise terms are independent and identically distributed. Since 
we require the trace of the noise covariance matrix to be equal in all the schemes, the covariance matrix for this 
scheme will be 



K 



iid 



Tr(K) 
N 



I. 



We denote this model as Scheme-iid indicating that the relay noise terms are i.i.d. Again, the optimal relay 
amplification vector for this model can be obtained from Theorem 1 by replacing K with Thus we 

have 



diid = 1^ 



n -1 



(f^OgY- 



(21) 



The corresponding SNR is given by 

SNRiid = ^'^'i?(f0g)^ 



(GGtIf^P. + fftp + -^ll0l 



Tr(K). 



(fog) 



(22) 



2) Relays with no correlation knowledge: Scheme-10: In this scheme, correlation is induced between the relay 
noise terms. However the marginals remain the same as in Scheme-00, n^^ ^ CA/^(0, Kjj); i.e. noise at relay Ri 
is AWGN with variance Kjj. Note that Kjj represents the correlation between the noise terms ur. and nR. for 
i 7^ j, and in general may not be equal to zero. Since the relays do not utilize the correlation structure, the relay 
amphfication vector is the same as in the case of Scheme-00, i.e.. 



dio = doo- 



The SNR achieved with the relay gain Dio = diag(dio) is given by 



SNRio 



|(f0g)tD-i(f 0g)|VPK 

(f0g)tDro^ADro^(f0g) • 



(23) 
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As there is no correlation knowledge at the relays, the relay operation here involves only co-phasing of the input 
signal and does not involve noise cancelation. It is straightforward to see that Scheme- 11 will perform better than 
Scheme- 10 as it exploits the correlation structure in designing the relay amplification factors. Rewriting the results 
from Theorem [T] the optimal relay amplification vector that exploits full correlation knowledge is 

dn = K([(GKGtpR + fFtp + K)]"'(f 0g)^ . (24) 

and the SNR achieved with dn is given by 

SNRii = PPR{i g)t (gKG^ Pr + fftp + kY^ (fog). (25) 



(f 0g)^(G(K0l)GtPR)-i(f 0g) 

SNRio = 1 (27) 

(f g)^(G(K I)GtP^)-i(GKGtPjj)(G(K I)GtP^)-i(f g) 



It is not clear whether Scheme- 10 is also inferior to Scheme-00 and Scheme-iid. It is also not clear which is 
better among Scheme-00 and Scheme-iid. Answering these questions will provide useful insights into the impact 
of correlation in multi-hop AF relay networks. We address this in the rest of this section. 

B. Asymptotic analysis 

We compare the schemes for the two extreme cases of Pr —>■ cxd and P oo. 

1) Relay Power Pr ^ oo: From the property of the optimal network at high Pr, we have from ([TOll 

SNRii = f^K-^fP. (26) 
At very high Pr, Dio = G(K I)G^ Pr and A = GKG^ Pr . Substituting these in we obtain 

Pj?)-i(GKGtPjj)(G(K0l 
(ft(K0l)-^f)^P 

ft(K0l)-iK(K0l)-if ^ ^ 

It is evident from (l20b that the SNR of Scheme- 11 at very high relay power is 

SNRoo = Pft(K0l)-if. (29) 

Similarly the SNR of Scheme-iid is readily obtained as 

ftfP 

SNRiiH = , s , ■ (30) 

Trace(K)/iV ^ ' 

Remark: Except for Scheme- 10, all the relay schemes turn out to be a point to point SIMO channel. It must be 
noted that the noise at the multiple antenna receiver is correlated and the noise covariance matrix is the same as 
the relay noise covariance matrix of the original network. For Scheme- 10, the multiple antenna receiver is unaware 
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of the noise correlation and assumes the noise covariance matrix to be K I instead of K. 

To answer the question whether correlation hurts when there exists no knowledge of it, let us consider the 
difference SNRio — SNRqo- For the simple case of = 2 and for real channels, we have 

-2/i/2Ki2(/?Kii + f^K22) 

SNRio - SNRoo = ' f2T^ lof fJ ^^l) 

/iKii + /|K22 + 2/1/2K12 

where f = [fi f2]^- The difference can be either positive or negative depending on the term — 2/1/2K12. For 
example, when the signals are positively correlated (/1/2 > 0) and the noise components are negatively correlated 
K12 < 0, then correlation helps as part of noise gets canceled. Similarly when both the signal and noise components 
are correlated in the same direction (positive or negative) then correlation hurts as it increases the noise power 
Therefore the overall effect of correlation can only be determined from the average behavior. 

Theorem 3: At high relay power Pr, Scheme-10 outperforms Scheme-00 in terms of both average SNR and 
average rate. That is 

E [SNRio] > E [SNRoo] and E[Rio] > E[Roo] 

where Rio = log(l + SNRio) and Roo = log(l + SNRoo). 
The proof is presented in the Appendix. 

The above result is significant as it suggests that, in the average sense, correlation does not hurt even if correlation 
knowledge is not available. To determine the relationship for the rest of the cases, we take the expectation of SNR 
in (126]|, dH, and 

N N 

E[SNRii] = P^E[|/i|2](l/Ai) = Pj^l/A. (32) 
1=1 1=1 

N N 

E[SNRoo] = P^E[|/i|2](l/K,i) =Pj^l/K,i (33) 

i=l i=l 

where E[|/j|2] = 1 while Aj and K a are the i*^ eigen-value and i^^ main-diagonal term of K respectively. E[SNRii] 
is greater than E [SNRoo] follows from the reason that the ordered vector containing the eigen-values of K majorizes 
the main diagonal of K, i.e. A ^ d. We also have E [SNRoo] > E[SNRjjjj] due to the harmonic-arithmetic mean 
inequality. Therefore we obtain the following relationship: 

E[SNRii] > E[SNRio] > E[SNRoo] > E[SNRjj^] (35) 

Thus in the high relay power regime, correlation does not hurt even when the system is ignorant of the underlying 
noise correlation structure. Importantly, it is also clear that there is a performance gain when the correlation structure 
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2.2 




Fig. 5. Average rate of a single user system as a function of total relay power Pr for P = Pi = 10. 



is exploited. 

2) Source Power P oo: At very high source power, the optimal relay network (Scheme- 11) is equivalent to 
a MISO system. From (fTTI ). we have 

SNRn = PRg^g. (36) 
Similarly, it can be shown that the rest of the schemes also achieve the SNR in (l36l ) at high P. 

SNRio = SNRoo = SNRjid = P/jgtg. (37) 

It can be noticed that this scenario represents the maximal ratio transmission (MRT) where the relays act as a 
multiple antenna transmitter. This is intuitive as the noise at the relays is negligible compared to the source power 
and therefore the effect of the noise covariance matrix is non-existent. We now verify the results with numerical 
analysis in the following subsection. 

C. Numerical Results 

We consider a common-interference based model to generate the relay noise covariance matrix. Each relay in 
addition to its thermal noise observes faded versions of common interference signals. The effective relay noise 
including both the interference and the local noise is given by 

Q 

Heff = ^ hkifc + Hr (38) 
k=i 
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10 20 30 

Source Power P in dB 



Fig. 6. Average rate of a single user system as a function of source power P for Pr = Pi = 10. 



where is the signal transmitted by the k^^ interferer and hk is the channel between the relays and the fc*^ 
interferer. The covariance matrix is given by 



Q 



K = E[neff n;^ff] = hkh^P/, + I. 



(39) 



k=l 



The total interference power is Pi = Ylk=i ^ik- Fig- [3 plots the average rate per channel-use as a function of 
the total relay power for P = P/ = 10 for the case of two relays and one interferer. The main observations are 

1) The performance order of the schemes in (l35l) which was obtained at high Pr is valid at all values of Pr. 

2) The difference {Ru — i?io) which indicates the benefits of learning correlation increases with relay power. 
With increase in relay power, the schemes diverge in performance, which is an indication that correlation 
impacts more at high relay power. 

In Fig.|6l we plot the average rate as a function of source transmit power for Pr = Pj = 10. Clearly, with increase 
in source power, the schemes achieve the MRT bound. Therefore the effect of correlation is less pronounced with 
increasing source power. Fig. |7] shows the average rate as a function of the interference power for one interfering 
node and two relays. As one can expect, the average rate decreases with P/. However Ru does not reach zero 
even at infinite interference power. This can be explained through the following: With one interfering source and 
N relays, the average SNR for Scheme- 11 at very high Pr is given by 

1 



E[SNR 



iij 



(iV- 1) +E 



(40) 



_1 + P/ II hi ||2_ 

This suggests that the average SNR is at least — 1 irrespective of the interference power. However, the average 
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Total Interference Power P in dB 



Fig. 7. Average rate of a single user system as a function of interference power Pi for P = 10 and Pr = 100. 

SNR decreases when the number of interfering sources increases. This is due to loss in degrees of freedom due to 
increase in the range space of interference. 

Fig. [8] shows the average sum rate versus the number of relays. The covariance matrix is generated with the help 
of 9 interferers with total power Pj = 200. It can be noticed that the sum rate increases at a greater rate when the 
number of relays is greater than the number of interferers. As Pj oo, Ru when the number of relays is 
less than the number of interferers {N < Q). However for (A^ > Q), Ru does not vanish even when Pj oo. 

Remark: As we discussed earlier, the relay design problem is a tradeoff between canceling the interference and 
maximizing the signal power. At very high interference power, it is important to cancel the interference. At high 
Pr since the network behaves as a SIMO system with antennas at the receiver, up to {N — 1) interfering sources 
can be rejected. When the number of interferers is more than the total number of relays, interference cannot be 
completely nulled. 

D. Impact of channel knowledge at the relays 

Throughout this paper, we assume that the relays have perfect knowledge of all the channels in the network. It 
may seem that complete channel knowledge is required at the relays for the relationship in ( [35l ) to hold. In Scheme- 
11, since the relays know the noise covariance matrix, they must also know the channels between the relays and 
interferers (when the number of interferers is less than the number of relays). This information is not available at 
the relays for Scheme- 10. Still, it performs better than Scheme-00 where there is no noise correlation. In fact, it 
can be shown that Scheme- 10 outperforms Scheme-00 even when there is absolutely no channel knowledge at the 
relays including the source-relay and relay-destination channels. This is also true when the relays only have local 
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Fig. 8. Average rate of a single user system as a function of the number of relays with 9 interferers of total power Pi = 200 (23 
dB),andP = 10. 

channel information. That is, relay Ri knows only its forward and backward channels and g(z). The following 
theorem states these results: 

Theorem 4: The average achievable rate for a relay network is lesser with relay noise covariance matrix K 1 
than with any general K, where only local channel knowledge is available at the relays. In other words 



Further, when there is absolutely no CSI available at the relays, Scheme-10 outperforms Scheme-00 in terms of 
the average SNR, i.e. E[SNRio] > E[SNRoo]. 



It is now clear that relay noise correlation is always helpful regardless of the channel and correlation knowledge 
at the relays. For maximizing signal power channel knowledge is essential while correlation knowledge is required 
to minimize interference power. With perfect channel knowledge, increasing the number of relays (while keeping 
the total relay power fixed) is helpful as signal power is increased due to coherent combining. 



In this section, we consider another application of the relay network with correlated noise. Consider a three- 
hop relay network as shown in Fig. |9] This is an example of a multi-stage relay network. There are N and M 
relays in the first and second stages respectively. Let the power expended by the source, the first stage relays, 
the second stage relays be Pq, Pi and P2 respectively. We denote this network with the shorthand notation 




(41) 



Proof: Refer to the Appendix. 



IV. Three-hop parallel relay network 
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Fig. 9. Three-hop parallel relay network 



P2 



9^3(^0) di, Pi, H, d2, P2, g)- In the first slot, the received signal at the first stage relays is given by 



ri = fx + ni. 



(42) 



The relay transmitted signal is given by = Diri, where Di satisfies the power constraint of the first stage 
relays, i.e. Tr(DiD|(fftp + I)) = Pi. Here Di = diag(di). In the second time slot, the relays at the second stage 
receive 

ra = HDiri + n2 = HDifx + HDini + n2. (43) 



In the third slot, the signal transmitted by the second stage relays is t^^ = D2r2- The power constraint at the 
second stage relays is given by E [H = P2, which results in 



Tr (D2D^(HDifftDlHtPo + HDiDlnt + I)) = P2. 



(44) 



Notice that the choice of Di affects the power constraint at D2. The received signal at the destination is given by 



y = g^tr, = g^D2(HDi(fx + ni) + n2) + n 

= g^D2HDifx + g^D2HDini +gTD2n2 + n. 



(45) 
(46) 



Note that both the relay stages operate at their maximum power. For any feasible D2 and Di, the SNR is defined 
as 



SNR(Di,D2) 



Ig^DaHDifrPo 



gTDsHDi ||2+ II gTD2 ||2 +1' 



(47) 
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The relay design problem is stated below as 



max SNR(Di,D2), (48) 

D1,D2 

s.t. Tr(DiDl(fftp + I)) = Pi, (49) 

and Tr (D2D^(HDiffTDlH^Po + HDiDIh^ + I)) = Pa- (50) 



To solve this problem directly can be difficult as the constraints are interdependent. We therefore approach 
the problem differently. Now consider any feasible Di for the first stage relays. Once the relay operation for 
the first stage relays is fixed, the network reduces to a two hop parallel relay network with correlated noise 
^2{Po, fnew , D2 , ^^2 , g) where 

fnew = HDif. (51) 
By comparing (1461) and ([2]), relay noise in the reduced two hop network can be found as 

nR = HDini+n2. (52) 
Note that the noise in (|52l ) is correlated and the corresponding covariance matrix is given by 

K = HDiDIh^ + I. (53) 

Now that the network is a dual-hop parallel relay network, we can utilize the result in Theorem [T] to solve the relay 
optimization for the second stage relays. We have from ([7]) 

d^(Di) = K(^[GKGtP2 + (fnew4w0 1)^0 + K 01]"' (fnew 0g)) • (54) 

We therefore are able to determine the optimal design for the second stage relays for all feasible relay gain 
vectors for the first stage. Similarly if we can find optimal Di for all possible D2 we can then resort to iterating 
between the solutions. However the problem is that the second stage relays' power constraint is affected by Di. 
Therefore fixing D2 and finding the optimal Di can violate the power constraint for the second stage. One way 
to solve this problem is to use the reciprocity property of AF relay networks recently proved in [2]. We state the 
theorem from [2] in the following. 

Theorem 5: [2] The capacity of a three-hop relay network is unchanged when the role of the transmitter and 
the receiver is switched while maintaining the same transmit power at each hop for both the original and reciprocal 
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TABLE I 

Algorithm: three-hop Relay Network Optimization 

1) Choose any Di that satisfies the first stage power constraint. 

2) Reduce the network to a two hop network with noise correlation. 

3) Find the optimal D^(Di). 

4) Begin Loop: 

5) Determine the reciprocal network. 

6) Optimize the second stage relays. 

7) Repeat until convergence. 

8) End Loop 



channels. This is equivalent to 

9^3(^0, f,di, Pi, H, d2, ^2, g) =9^3(i^2,g,'^2d2, Pi, HT,Atidi,Po,f), (55) 

where the constants ki and K2 ensures power constraint compliance in the reciprocal network. (Note that we allow 
some abuse of notation as the above statement implies that the capacity of the two networks are the same.) 

Using the above result, one can reverse the direction of the communication in the relay network. For the reciprocal 
network, the gain of the first stage relays is related to the gain of the second stage relays through a constant 
multiplying factor which can be obtained from the power constraint as follows: 

9 Pi 

= — r (56) 

Tr D2DT(ggtP2 + I) 



Similarly ki is found as 



2 ^^rjs 

Ki = — 7 ; ; ^. (57) 



Tr (DiDi(K2HTD2ggtD^H*P2 + k2htD2D^H* + I) 



As discussed earlier, the reciprocal network can be reduced to a two hop network and the relay gain for the second 
stage can be optimized. Using Theorem[T] we can find the optimal gain for the second stage of the reciprocal network 
(which is proportional to the first stage relay of the original network) for every feasible K2D2. The procedure then 
will be to find the optimal relay gain for the second stage having the first stage gain fixed. The next step is to find 
the reciprocal network and optimize the second stage, and iterate this procedure till convergence is achieved. The 
algorithm is listed in Table Jl 

Remark: Notice that the capacity of the network increases with every iteration as each iteration involves 
optimization. This suggests that the algorithm will eventually converge. Although it is not known if the algorithm 
guarantees a global optimum solution, based on extensive numerical simulations, we conjecture that the algorithm 
outputs the global optima whenever there is convergence. 
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Fig. 10. Pictorial representation of the relay optimization algorithm. 

Fig. [TT] shows the convergence of the algorithm for a variety of initial relay amplification factors. In the simulation, 
the channel matrices are 

f=[16], g=[4 -3], H=[2 -3;4 2]. 

We use the following initial inputs for the first stage relay gains every time the algorithm begins: 

di(l) = [1 0], di(2) = [0 1], di(3) = [-2 1], di(4) = [2 1], di(5) = [-20, - 1]. 

The transmit power at each hop is unity (Pq = Pi = P2 = 1). At the start of the algorithm, the inputs will be 
scaled to satisfy the power constraint. For each input, the progress of the algorithm (measured in terms of SNR at 
the destination) is displayed in Fig. [TT] For all these inputs, the algorithm converges to an SNR value of 6.5638 
while the first and second stage relay gains converge to 



-0.0823 - 0.1633], d'a = [0.2533 0.2566]. 
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Fig. 11. Convergence of the algorithm for different initial inputs (di). 

In the following we study the characteristics of three-hop relay networks in the high-power regimes. The proofs 
are available in the Appendix. 

A. Properties of 9^3 (^o , f , di , Pi , H, da , P2 , g) 

1) P2 —>■ cxd: When there is no power constraint (or infinite power) for the last stage relays, the three-hop network 
for the optimal D2 reduces to a two hop network with the destination consisting of multiple antennas. This 
is equivalent to 

$n3(Po,f,di,Pi,H,d2 = d^,P2 = oo,g) = 9l2(Po,f,di,Pi,H). (58) 

It is worth noting that the noise at the multiple antenna destination of the reduced network is correlated. 

2) Pq ^ 00: For the case where the source power is infinite, it is clear from the properties in Section Hl-CI that 
the noise terms at the first stage relays are negligible. This reduces the network to a two-hop relay network 
with a multiple antenna source. 

9^3(Po = 00, f , di = d°i, Pi, H, d2, P2, g) = ^2(A, H, d2, P2, g) 

The network optimization for the reduced network involves joint precoder (transmit covariance matrix) and 
relay design which are interdependent. We can utilize the algorithm in Table U for this optimization. 

3) Po,P2 —>■ 00: When the power of both the source and the last stage relays tend to cxd, the relay network 
reduces to a point to point MIMO channel with the constraint that the rank of the input covariance matrix is 
one. 

$n3 (i^o = 00 , f , di = d°i , Pi , H, d2 = d^ , P2 , g) = % (Pi , H) 
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Here £Hi (Pi , H) indicates that the network is a single-hop point to point channel. Thus the three-hop network 
optimization problem is equivalent to SNR maximization problem in the point to point MIMO channel. The 
optimal strategy is beamforming along the principal eigen-vector of the channel matrix. 

4) Pi, P2 00: For this scenario the last two stages can be removed and the network reduces to a point to 
point SIMO channel with receive combining vector §"""021101. Optimizing one set of relays is enough to 
achieve the point to point channel capacity. 

$K3(Po, f , di, Pi, H, d2 = d^, P2 = 00, g) = ^i(Po, f). 

0li(Po, f) indicates that the network is a single-hop point to point channel. 

5) Pq, Pi — >^ 00: When the source and the first stage relays have very high power, the network reduces to a point 
to point MISO system with transmit precoding vector D2HDif. Like the previous case, optimizing one set 
of relays suffices to optimize the whole network. 

$K3(Po = 00, f , di, Pi = 00, H, d2 = d^, P2, g) = ^l(P2, g) 

Remark: In all the above cases, if the transmit power at any hop is very high, noise at the receive side is 
negligible. When this occurs at the first or the last hop, (for example, Pq or P2 00) the hop can be removed. 
Here any rate that can be achieved in the original network is achievable in the reduced network and vice versa. For 
example consider the three-hop case where the source has one antenna, and 5 relays in the first stage. When the 
source power is very high (Pq 00) then the first hop network can be removed, and the reduced network consists 
of two hops. In the new network the source consists of 5 antennas and transmits a one dimensional signal. It is 
interesting that the two networks are equivalent. 

V. Conclusion 

In this work, we considered an AF relay network wherein the relay noises are correlated which may be due to 
common interference or multi-hop AF relaying. We obtained closed form expressions for optimal rate maximizing 
relay gains and maximum achievable rate when correlation knowledge is available at the relays for both single-user 
and multi-source scenarios. Further we showed that correlation does not hurt irrespective of channel and correlation 
knowledge at the relays. We also showed that correlation knowledge results in significant performance improvement. 
Analytical and simulation results demonstrate significant rate enhancement when correlation knowledge is exploited. 
With appropriate gains, the relays can perform distributed interference cancelation when the number of relays is 
greater than the the number of interferers. The performance improvement increases with the total relay power and 
the number of relays. As there are significant benefits in learning correlation, practical schemes to communicate 
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the correlation structure to the relays need to be explored. 

For three -hop AF relay networks, we addressed the relay optimization problem. By fixing relay gains for one of 
the stages we reduced the network to a two hop network with noise correlation. Combining the result we obtained 
for two hop networks with correlated noise and a multi-hop duality result from [2], we propose an iterative algorithm 
that increases capacity for every iteration. This algorithm can be utilized for joint optimization of relay gains and 
precoder in multi-antenna sources. 

Appendix 

A. Proof of Theorem |7] 

From the destination received signal in Q, the achievable rate is given hy R = log2(l + SNR) where 

dTGfftGtd*P 

It is worth noting that d should satisfy the sum power constraint of the relays. 

(fftp + k) l] d* = Pr (60) 

Now we are interested in maximizing the sum rate over all possible relay amplification vectors that satisfy the relay 
power constraint. The above problem is equivalent to maximizing SNR with respect to d. 

SNR° = max SNR 
d 

Incorporating ( [60l ) into the denominator of (|59l ). we have 

d^Gff^G^d*PPR 



SNR 



Define the matrices A and B: 



dTGKGtd*Pfl + dT [(fTtP + K) I] d* 

dTGfrtGtd*ppR 

dT [GKGtPR + (fftP + K) I] d* 



A = GKG^Pr + ( fftp + K ) 1 



(61) 



(62) 



B = GfTtG^PPn (63) 



Notice that the matrix A is Hermitian and positive definite while B is Hermitian and positive semi-definite. By 
Cholesky decomposition A = LL^, we have 
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Let w"*^ = d^L, we have 



SNR 



w 



Tl iBL tw* 



(65) 



The relay optimization problem is therefore 



SNR° = max 



w 



TL-iBL-tw* 



By utilizing Rayleigh Quotient we have 



SNR° = A^ax(L-iBL-t) 



(66) 



(67) 



w° = fv^ax(L-^BL-t) 



(68) 



where Amax(P) and i;max(P) are the principal eigen-value and eigen-vector of the matrix P. 
Since B = (f g)(f g)^PPR, we have 

SNR° = PPRAn,ax(L-^(f 0g)(f 0g)^L-t) 

= PPflA^ax((f 0g)^L-tL-i(f 0g)) 
= PPfiA^ax((f 0g)^A-i(f 0g)) 
= PPr{{ Qg)^A-\fQg) 



(69) 
(70) 
(71) 
(72) 



where dTOl ) follows from the property that the non-zero eigen values of AB are equal to those of BA. dTTI ) holds 
since L^^L^-*^ = (LL^)^-*^ = A^-*^. The optimal relay amplification vector is given by 



(73) 
(74) 
(75) 



d° = (L-tv^,,(L-iBL-t- 
= K(L-tL-i(f0g) 
= K(A-i(f0g))*. 

where k ensures compliance with the relay power constraint. Hence proved. 



B. Proof of Properties in Section I7/-CI 

1) Pr ^ oo: At high relay power Pr, for any relay amplification vector d, the destination noise is negligible. 
In that case, the input-output relationship can be expressed as 



y = (d g)'^fx + (d g)'^nR, 



25 

which is equivalent to a SIMO system, y = fx + nR , with the receive combining vector d0g. Now consider 
the maximum achievable SNR (corresponding to d°) from which is 

SNR° = PPij(f0g)tA^i(f0g). 

At high Pr, A = (GKGt)PR which leads to 

SNR° = P(f g)t(GKGt)-^(f g) = Pf^K-if . 

At high Pji, the optimal relay functionality becomes 

d° = K(GKGl')~^(f g) = KG^^K-^G^^Gf = KG^^K^^f. 

Thus d = KG^^K^-'^f achieves the MRC bound when Pr —>■ oo. 
2) P — > cxd: When the source power is high, it can be observed that the relay noise terms are negligible. The 
input-output relationship is then given by 

T (d f ) 



which represents a MISO system with the beamforming vector ^f^J^fy j^^Qf^ • Since (d f)^(d f)P = Pr, 
the input transmit power is Pr. For the optimal relay design d° we have A = (ff ^ I)P, which results in 

SNR° = Pnif g)t(frt I)"i(f g) = P^gtg. 

The system thus behaves as a multiple antenna transmitter with a single antenna receiver. ■ 

C. Proof of Theorem |2] 

From the destination received signal in ( [T3l ). we can notice that due to relay amplification, the relay network is 
equivalent to a scalar MAC. The sum rate is given by 

R = log(l + SNR) 

where 

d^G (ELifkfilP*) G*d' 
^''''= dTGKGd- + 1 ■ <'« 



Notice that d satisfies the sum power constraint of the relays. 



^fk^Pfe + K 01 



Pr (77) 
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Now we are interested in maximizing the sum rate over all possible relay amplification vectors that satisfy the relay 
power constraint. The above problem is equivalent to maximizing SNR with respect to d. 

SNR° = max SNR 
d 

Incorporating (TTTI ) into the denominator of (TtSI ). we have 

SNR= ^ — 



dTGKGtd*PR + dT 



Etifkfi^^fe + K)0l 



d* 



d^G (eLi fkfi^^fc) Gtd*P^ 



dT 



Define the matrices A and B: 



GKGtP^ + (J:;:^^ fkfTPfc + Kj I 

GKGtPn + (^Yl fkfk^fc + I 
B = G I^J^fk^Pfc^ G^Pr 



d* 



Then, proceeding as in the proof for Theorem 1, we have. 



SNR° = A^ax(A-iB) 



and the optimal relay amplification vector is given by 



d° = (v^ax(A''B))* 



Hence proved. 

D. Proof of Theorem \3\ 

From (l28]l and (l29]l we have 

where 



SNRoo = a SNR 



10 



a + b 



N 

1=1 
N N 



m 

K, 



i=l j=l j^i 



(78) 



(79) 



(80) 



(81) 



(82) 



(83) 



(84) 



(85) 
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Now consider the the difference in SNR of Scheme- 10 and Scheme-00, which is given by 

SNRio - SNRoo = ^ (86) 
a + b 

Notice that the difference SNRio — SNRqo is convex in the correlation term b as its second derivative w.r.t. b is 
non-negative. We assume that the magnitude and phase of the components of f are independent where the magnitude 
follows Rayleigh distribution while the phase is uniformly distributed in the interval [0, 27r]. Conditioned on the 
magnitude of the channels, the average of b over all channel phases is 0. By applying Jensen's inequality we have 

E[SNRio - SNRoo] > 0. 

Let us now consider the difference in rate Rio — Rqo- 

Rio — Rqo = IE 



, , 1 + SNRio 
log 



1 + SNR, 



00 



log lif . (87) 



It can be shown that the difference i?io — -Rqo is convex in b. Therefore through Jensen's inequality 

E[i?io - Roo] > 0. (88) 
Hence proved. ■ 

E. Proof of Theorem |?] 

For any constant relay gain D for the relays, the SNR of the two schemes where the relay noise covariance 
matrix is K and K I respectively are 



SNR 



10 



SNR, 



|g^Df|V 
gTDKD*g* + 1 



00 



gTD(K0l)D*g* + 1 

Note that the power expended by the relays is the same in both the schemes. When the relays have only local 
channel knowledge, the phase of the relays should be such that the phase of the forward and backward channels 
are canceled. We therefore have 

SNRio = and SNRqo = t 

b + a b 

where 

N N 

« = E E g«g*(j)d(i)d*(j)K., 

i=l j=l j^i 
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N 



b = J2\m'\d{i)\'K,. 



i=l 



c= (^^|d(i)g(i)f(i)|j 



The term a depends on both magnitude and phase of the channels while b and c depend only on the magnitude 
of the channels. It can be shown that log(l + SNRio) is convex in a. Given the magnitude of the channels, the 
average of a over channel phases is zero. Therefore by Jensen's inequality 

Ea,b,c [log(l + SNRio)] > Ef,,c [log (l + ^)] = Efe,c [log(l + SNRqo)] • 

Now consider the case where the relays have no CSI. Here the relay gain D is independent of f and g. Taking 
expectation of SNRio and SNRqo over f, we have 

II S^D PP 

where ¥,[f£^ = I. Observe that ( [89l ) is dependent on the phase of g while ( [90l ) is independent of phase of g. From 
Jensen's inequality it can be easily seen that E[Ef [SNRio]] > IE[Ef [SNRqo]]- Hence proved. ■ 

F. Proof of Properties in Section l/V-AI 

From (l46l ). the received signal at the destination for the three-hop case is given by 

y = g'^DaHDifx + g'^DaHDini + g^Dsna + n. (91) 

Let Di = y/PiD^ and D2 = y/^Y)2 where DH expends unit power at the i^^ stage relays, i G {1,2}. Similarly 
let X = \/7ox' where E[|x'p] = 1. Now the effective received signal at the destination is 

= g-D^HD^fx' + + + - . (92) 



The modified power constraints are 



Tr 



Tr(^D^D^t(^fft + -l^^ =1, (93) 
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1) Pa ^ oo: 

At very high P2, the normalized input-output relation becomes 



= g^D^ f HD^fx' + + ^ V (96) 

V vPo vPi/ 

This represents a two-hop hop network with multiple antennas at the destination where g^D2 is used as its 
receive combining vector. It is worth noting that maximal ratio combining (MRC) at the destination of the 
reduced network is capacity optimal. Let wmrc be the optimal receive combining vector. Now, we need to 
show that there exists a diagonal matrix in the original three-hop network such that g^D2 = cwmrc 
for some constant c. This can be achieved when D2(ii) = ^■^^^^^^^^ where k satisfies the power constraint 
of the second stage relays given in (|94l ). (Here D2(«^) denotes the i^^ element of the diagonal matrix D2.) 
Therefore we have 

fH3(^o,f,di,Pi,H,d2 = d°,P2 = oo,g) =$n2(Po,f,di,Pi,H). (97) 

2) Po ^ 00: 

We prove this by using the previous property and the reciprocity property of AF relay networks. 

^R3(Po = 00, f , di = d°i, Pi, H, d2, P2, g) = 9^3 (P2, g, ^J2d2, Pi, H^, AJidi = Kid°i, Pq = 00, f) (98) 

= ^H2(P2,g,K2d2,Pi,HT) (99) 

= lH2(Pl,H,d2,P2,g) (100) 



where ( 1981 ) follows from the reciprocity of multi-hop AF networks. ( |99l l follows from the preceding discussion 
on P2 00. (llOOl ) uses the reciprocity property again. 

3) Po,P2 ^00: 

For this scenario, we can neglect ni and n from ( [92l ) which results in 

y' = g^D^HDifx' + (101) 

V Pi 



g^D^ ( HDJ'fx' + ) (102) 



^P 



with the power constraints being 



II D^f f= 1 (103) 
Tr (b^B^^ (uD^S^n^ni + -1^ ^ = 1. (104) 
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( 11021) represents a point to point MIMO channel with D"f and g'''^D2 as its transmit precoding and receive 
combining vectors respectively. Notice from ( |103b that the precoding vector D"f has unit norm. Therefore 
the transmit power for the channel in (11021 ) is Pi. Now let us consider the point to point MIMO channel 
Y = HX + N. For any precoding vector u and receive combining vector v in this channel, there is an 
equivalent relay gain for the first stage and second stage relays in the original network such that the two 
systems are equivalent. That is, we can find D" and D2 such that the following conditions are satisfied for 
some constant c: 

D^f = u and g'^D^ = cv 

Since the precoding vector u has to be unit norm, the power constraint for the first stage relays is automatically 
taken care of, while the constant c takes care of the power constraint for the second stage relays. Therefore 
the original three-hop relay network can be reduced to a single hop point to point MIMO channel as given in 
the following. 

^3(^0 = 00, f , di = d°, Pi, H, d2 = d°, P2 = 00, g) = lHi(Pi, H) 

4) Pi,P2^oo: 

At very high Pi and P2, the normalized input-output relation for the three-hop network becomes 

= g^D^HDifx- + g^D2iroyni ^^^^^ 

vPq 

= gTD^HD^(fx' + -^). (106) 

The power constraints are 

Tr(^DMt(^fFt + -l^^ =1, (107) 

Tr ^D^D^^ ^HD^fftD^H^ + HDi^iHt ^ ^ ^ ^ ^^^^^ 

(1106b represents a point to point SIMO system y = fx + n with transmit power Pq and g^D2HD" as its 
receive combining vector. We now have to prove that there exists diagonal matrices D" and D2 satisfying 
the power constraints such that g^D2HD" = cv for some non-zero constant c. Here v is the optimal MRC 
vector for the point to point SIMO system. This can be achieved through the following: First, fix any D2. 
Then, find D" such that g^D2HD" = cv where c ensures compliance of (1107b . Now scale D2 to meet the 
constraint in (1108b . Hence the original three-hop network is equivalent to a point to point SIMO channel and 
is represented as 

lR3(i^o,f,di,Pi,H,d2 = d^,P2 = oo,g) = 9^i(Po,f)- 
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5) Po,Pi^oo: 

By reciprocity, the proof for the previous case holds for this case as well. ■ 
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